Holstein light quantum polarons on the one- dimensional lattice 



O. S. BarisicE 

Institute of Physics, Bijenicka c. 46, HR- 10000 Zagreb, Croatia 

The polaron formation is investigated in the intermediate regime of the Holstein model by using 
an exact diagonalization technique for the one-dimensional infinite lattice. The numerical results 
for the electron and phonon propagators are compared with the nonadiabatic weak- and strong- 
coupling perturbation theories, as well as with the harmonic adiabatic approximation. A qualitative 
explanation of the crossover regime between the self-trapped and free-particle-like behaviors, not 
well-understood previously, is proposed. It is shown that a fine balance of nonadiabatic and adia- 
batic contributions determines the motion of small polarons, making them light. A comprehensive 
analysis of spatially and temporally resolved low-frequency lattice correlations that characterize 
the translationally invariant polaron states is derived. Various behaviors of the polaronic deforma- 
tion field, ranging from classical adiabatic for strong couplings to quantum nonadiabatic for weak 
couplings, are discussed. 



I. INTRODUCTION 



Charge transport observed in the limit of the low 
electronic density is characterized by significant pola- 
ronic effects when the electron-phonon coupling is strong 
enough. Indeed, the effects of this kind have been 
reported recently in a number of experiments on sys- 
tems ranging from colossal magnetoresistive manganites 
to high-T c superconductorsji*^ which are currently the 
subject of great scientific interest. In this context, a bet- 
ter understanding of the polaron states is of particular 
importance. The local polaronic correlations and the po- 
laron dispersion properties are characterized by adiabatic 
and nonadiabatic contributions combined with the effects 
due to the discreteness of the lattice. In this respect 
the polaron problem contains most of the conceptual ele- 
ments important for the dynamics of other problems with 
localized solutions. In particular, the theory of quantum 
polarons has its analogies in the more general context of 
soliton physics, invoked for broad classes of materials. 

The experimental evidence of strongly coupled pola- 
ronic carriers is frequently provided in terms of the ther- 
mally activated d.c. conductivity, which is accompa- 
nied by a broad contribution to the mid-infrared opti- 
cal conductivity^ On the other hand, in photoemis- 
sion spectroscopy measurements, one observes a small 
electron quasi-particle weight at the Fermi level, with 
most of the spectral weight transferred towards high 
frequencies In order to rationalize such strong pola- 
ronic correlations, assuming a local electron-phonon in- 
teraction, one may apply the Holstein model. Within 
the Holstein model^ in the strong-coupling limit the 
electron-phonon spatial correlations reduce to a single 
lattice site and one obtains heavy, strongly-pinned, small 
polaron states characterized by a large lattice distor- 
tion. For such states the polaron binding energy e p 
determines the thermal activation energy, Ui c . « s p /2, 
and the optical activation energy ui op t ~ 2e p i^ Yet, even 
for quite strong couplings the adiabatic spreading of the 
lattice distortion over a few lattice sites might signifi- 
cantly reduce the pinning barrier for the polaron motion, 



^d.c. <C u>opt- In particular, in the one-dimensional Hol- 
stein model weakly-pinned small polarons with a consid- 
erable binding energy are obtained between the strong- 
and weak-coupling limits for moderate values of the adia- 
batic parameter. One of the key objectives of this work is 
to achieve a qualitative understanding of this interesting 
parameter regime, which can be relevant for unconven- 
tional small-polaron behaviors reported for various ma- 
terials. 

The accurate treatment of adiabatic correlations is a 
difficult problem for itinerant polarons. Several numer- 
ical methods have been recently suggested for the one- 
dimensional Holstein model in order to spectrally resolve 
the polaronic correlations in the single-electron and op- 
tical conductivity spectra i 1 ^ 11 ' 12 These studies clearly 
show the fundamental difference between the nonadia- 
batic and adiabatic regimes, although the reliability of 
the methods at low frequencies is not entirely clear. In 
the adiabatic regime one should obtain the excited coher- 
ent polaron bands with finite spectral weight below the 
phonon threshold, each of the bands corresponding to one 
of the adiabatically softened local phonon modes 42 Yet, 
the spectral weight assigned to such low-frequency adia- 
batic correlations seems to be reported in the literature 
only for the case of Diagrammatic Quantum Monte Carlo 
calculations pi^ The Dynamic mean-field theory (DMFT) 
applied to the Holstein polaron problem^ exact in the 
infinite dimensional limit, gives only one excited co- 
herent polaron band below the phonon threshold. For 
low-dimensional systems the approache D 1 ^ 1 ? like DMFT 
that treat the electron self-energy as a local quantity, 
£(fc,u>) — > provide approximate results. Namely, 

in order to describe correctly the local adiabatic correla- 
tions within the diagrammatic approach, the vertex cor- 
rections involving phonons at different lattice sites i ^ j, 
which contribute to the fc-dependent (nonlocal) part of 
the electron self-energy E(fc, u>), have to be taken into ac- 
count. In fact, the important vertex corrections turn out 
to be those within the range of the adiabatic correlation 
length^ This length is restricted by the distance over 
which the electron can move before the lattice distortion 
evolves significantly in time. 
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In order to obtain the possibly most accurate picture of 
the low-frequency polaronic correlations a previously de- 
veloped approac h 1 ^ 1 ? of exact diagonalization of the infi- 
nite lattice polaron problem is employed here. By provid- 
ing eigenenergies and corresponding wave functions, this 
approach, unlike any other developed up to now, gives the 
complete information about the coherent part of the po- 
laron spectrum below the phonon threshold. In addition 
to the polaron spectra and electron spectral properties, 
the phonon spectra that describe the displacement and 
momentum correlations of the lattice can be analyzed in 
this way. In fact, this work provides, to our knowledge, 
the first comprehensive study of lattice properties within 
the translationally invariant polaron theory. According 
to the regime of parameters, qualitatively different be- 
haviors of the lattice deformation field are obtained in 
the low-frequency range. Such findings can be of exper- 
imental interest, for example, in the context of recent 
neutron diffuse scattering (Huang scattering) or pulsed 
neutron-diffraction measurements^ 



II. GENERAL 

The Holstein polaron arises from the local interaction 
of the electron in the nearest-neighbor tight-binding band 
with the dispersionless branch of optical phonons. For 
the one-dimensional lattice the Hamiltonian reads 

H = -2i^cos(/c) c\c k +u ^b\b q 

k q 

- ff/\^5>U c *( 6 -9 +ft «)' (i) 

k.q 

where N is the number of lattice sites (N — > oo). c\ (cfc) 
and b J (b q ) are the creation (annihilation) operators for 
the electron state with momentum k and phonon state 
with the momentum q, respectively. Besides t/ui and 
g/oJo, the two other parameters used here are the bind- 
ing energy of the polaron in the atomic limit (t = 0), 
£ v = 9 2 / u o, and the dimensioneless parameter A = e p /t. 
Notice that e p and A, unlike g, are adiabatic parameters 
in the sense that they do not depend of the lattice mass. 
The phonon energy luq is taken throughout this work as 
the unit of energy. 



G K (t) = -i(0\f Mi)4]|0>, 

with cjf = l/VNJ2j e tKj c^ and T the time ordering 
operator, is evaluated for the ground state of the sys- 
tem |0) in which no electrons are present. As there are 
no hole excitations for such a state, the spectral func- 
tion that corresponds to the Fourier transform of Gk if) , 
A K iu>) = — 7r _1 ImG\R:(w), describes only electrons above 
the polaron ground state energy Eq, If Axi^j) is decom- 
posed into two parts with respect to the phonon thresh- 
old, 

A K (w) = A^iEo < u < uj + E ) + A^uj >uj + E ) , 

the low-frequency coherent part of the spectrum can be 
expressed as 

UJ<Eq+UJo 

A<{w)= £ A® 8(u-E$), (2) 

i 

with 

4? = i<*£i4>i a , 

where \c' K ) represents the free electron state without 

phonons, E^ is the energy and ) is the wave func- 
tion of the translationally invariant polaron states, and K 
and i are used to denote the momentum and the number 
of the polaron band, respectively. At zero temperature, 
for energies below the phonon threshold uj < E + u , 
nonelastic scattering is forbidden by energy conserva- 
tion, which makes the imaginary part of the exact elec- 
tron (and phonon) self-energy arbitrarily small. For this 
reason, the incoherent contributions can be neglected in 
Eq. J3J, and A^(uj) may indeed be regarded as a signa- 
ture of the coherent polaron spectrum, but only below 
the phonon threshold. 

B. Lattice properties 

The phonon propagator used here (not to be confused 
with the displacement-displacement correlation function) 
is defined as 



A. Electron properties 

In the present work different dynamical correlation 
functions, which describe properties of the system in 
the dilute limit, are investigated at frequencies below 
the phonon threshold. This subsection serves to outline 
the relevant notation while examining some general as- 
pects of electron correlations. The next subsection, on 
the other hand, deals with lattice correlations. 

The electron propagator, here defined by 



D Q (t) = -i(* |r [M*H]|*o> . (3) 

where |^o) ls the ground state of the electron-phonon 
system. The phonon spectral function Fq(o;) = 
— 7r _1 Im_Dg(tj) is even in Q and generally asymmetric 
in uj. Fq(u>) may be used to spectrally resolve the boson 
commutation relation 

y>OC 

/ [F Q (u)-Fq{-w)]du = l. (4) 
Jo 
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The mean number of phonons with the wave vector Q 
present in the ground state is given by 

N% = (* |^6q|*o> = f F Q {u) . (5) 

J — oo 

The spectral function corresponding to the displacement- 
displacement correlation function, 

D* Q *(t) = -i(* \f [A Q (t)A- Q }\* ) , (6) 
with Aq = b_Q + 6q, is written as 

p q X (u) = Fq(lj) + Fq{-w) + 2H Q {u) , (7) 

where Hq(lu), and thus Pq x (lu), is an even function in Q 
and u). Hq(lu) represents the spectral function related to 
the correlation function 

Dg(t) = -i/2 (* |T + 6 + _ (*)6S,]|*o> ■ 

It is straightforward to express the spectral func- 
tions for the displacement-momentum and momentum- 
momentum correlation functions in terms of Fq(u>) and 
Hq(us). For the 3 = ground state with no phonons, one 
obtains Fq{uj < 0) = and Hq(uj) = 0. 

The properties of the general expressions for Fq(u>) 
and Hq(lu) can be further specified for the polaron prob- 
lem. Just as in the electron case, the polaronic corre- 
lations manifest themselves in the formation of coherent 
polaron bands below the phonon threshold. Accordingly, 
the low-frequency coherent part of the spectral function 
Fq(u>) can be expressed in terms of translationally in- 
variant polaron states, 

F Q (w) = F<(\u\<u ) + F>(\u\>u ) 

\u\<u) 

F<H = F ±.Q S(wtE^±Eo), (8) 

i 

with 

= K*2Hl*o>| a , f% = |<*« |M*o>| 2 ■ 

Similarly, the coherent part \uj\ < uj of Hq{lo) is given 
by 

|w|<w 

ff<H = J2 H Q t("TE$±E ) (9) 

i 

H® = l/2<*o|M*Q><*g ) l&-Ql*o> + h.c. . 

In the absence of electrons, the phonon propagator is 
purely local, having a trivial temporal dependence, i.e., 



Fq{lo) = 6(oj — Wo). The phonon spectral weight for 
\w\ < ujq in Eqs. (jBJ and JSJ is therefore entirely due to 
the coupling with the electron, Fq (u) , Hq (w) ~ 1/N. 

In the long-wave limit Q — > the phonon spectra in 
Eqs. (jHJ and © are smooth functions of Q, including 
Q = 0. Otherwise, one would find that eigenstates below 
the phonon threshold involve phonons at infinite distance 
from the electron. Such phonon excitations uncorrelated 
to the electron are forbidden by the energy constraint. In 
Eq. the phonon Q — mode couples only to the total 
number of electrons J2k c k c k> which is a constant of mo- 
tion. The homogenous Q = part of the Hamiltonian Q 
can therefore be separated outi^i Its eigenstates corre- 
spond to a displaced harmonic oscillator of frequency ujq. 
Assuming one electron in the system, for \u\ < ojo in 
Eqs. (jHJ and 10 it is straightforward to obtain 

F® Q=0 = F% =Q = ffgl = 6 ifi s p /Nu . (10) 

The reference case, which illustrates well the funda- 
mental difference between the role of the lowest i = 
and the excited i > polaron bands in the lattice cor- 
relation functions, is the limit of self-trapped small po- 
larons (e p ^> u>o, A 3> 1). Namely, in this limit the 
polaron spectrum below the phonon threshold exhibits 
a few well-separated, narrow polaron bands. The lowest 
(i — 0) band contribution to the lattice correlations de- 
scribes the heavy polaronic lattice distortion, while the 
excited (i > 0) bands correspond to the local harmonic 
modes softened below u>o, and carried along by the mov- 
ing polaron. The phonon spectral weight related to the 
lowest (t = 0) band involves many phonons, which for 
small polarons, are equally distributed over all Q's, 

F% « F% « < « N tot /N , (11) 

with N tot sa e p /oj S> 1 representing the mean total num- 
ber of phonons in the ground state, 

N tot = J2(*o\blb Q \* ) . (12) 
Q 

Ntot is large here because of the large lattice distortion 
that accompanies the formation of the self-trapped po- 
laron. For |w| < W^°\ where is the bandwidth of 
the lowest polaron band, one finds from Eqs. 11 Oil and 
(1 1 H that the local lattice distortion of the self-trapped 
small polaron behaves classically for any Q, 

Fi%'F W Q ^l/N, (13) 

i.e., as if the boson operators in Eq. 10 commute. Hence, 
Fq(oj) as an even function of w for ui small describes 
the classical nature of the small polaron self-trapping. 
This represents the generalization of well-known results 
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for large adiabatic polarons^ to the small-polaron case. 
Furthermore, the distortion moves along the lattice on a 
long time scale, which is determined by the bandwidth of 
the lowest polaron band W^°\ -C <^o- Such motion 
is characterized by a small lattice kinetic energy, i.e., 

where the left-hand side corresponds to the spectral func- 
tion of the momentum-momentum correlation function. 

The total contribution of the zth polaron band to the 
phonon spectral weight can be analyzed in terms of 

*f = ^ = £<- 

Q Q Q 

(14) 

Unlike for the lowest (i = 0) band, the contributions 
itTHl for the excited (i > 0) bands are independent of 

parameters; m 1 and F_ w H (l > w 0. The contri- 
bution of the excited bands appears only in Fq(u>) at the 
lo > side, which is characteristic for harmonic modes 
that are carried here by the moving polaron, again in 
generalization of the results obtained in the context of 
large adiabatic polarons^ to the small polaron case. 



C. ET method 

The polaron states below the phonon threshold, which 
define the low-frequency electron and lattice correlations 
discussed in last two subsections, can be calculated prac- 
tically exactly by applying the exact-translational (ET) 
method. This method, introduced in Ref. ITflL for small 
and moderate values of the adiabatic parameter t/u>o 
gives the lowest polaron band for the infinite lattice prob- 
lem with the highest accuracy^ In addition, the ET 
method as implemented in Ref. ^3 reveals, with similar 
accuracy, the excited polaron bands below the phonon 
threshold. It should be mentioned that the current work 
provides only a brief overview of the ET m ethod, where as 
detailed discussions can be found in Refs. Il. c llii24l24 

The ET method resolves the spectrum below the 
phonon threshold in terms of the basis set 

\n) K = -^=V]e lJS ' J0 ct|Tio,n_i ) Tii ) ... ) n m) ...)j , (15) 

V ./V . 
j 

where n m defines the number of phonons at a given dis- 
tance m from the electron. On the left hand side of 
Eq. ltl5)i n is used as an abbreviation for the whole set 
of quantum numbers n m . The translational symmetry of 
the problem JQl is exploited by the construction of the 
states l(T5|l . which are the eigenstates of the total mo- 
mentum. In calculations, the basis set l(T5|l is kept finite 
by restricting the number of phonons and their distance 



from the electron. Such truncation of the basis is justi- 
fied for the lowest part of the spectrum below the phonon 
threshold, which is characterized by the finite electron- 
phonon correlation length d. Indeed, in the weak cou- 
pling regime, d is nonadiabatic, independent of g, and 
determined by the ratio t/uj . For strong couplings d is 
independent of the lattice mass; d ~ A -1 . 

In the representation l(T5)l . the Hamiltonian is given 
by a complex Hermitian matrix. However, with the ap- 
propriate transformation of the basis J loll , similar to the 
one proposed for the antiferromagnetic chain problem in 
Ref. Hi 

\n+) K = {\n) K + \n) K )/V2 (16) 
\n-) K = {\n)K-\n) K )/iV2, (17) 

the Hamiltonian takes the form of a real symmetric ma- 
trix. With respect to the previous implementations, this 
represents a further improvement of the ET method due 
to the considerable reduction in computational effort as- 
sociated with it. In Eqs. ijlfijl and itTTjl \n)K is ob- 
tained from \n)K when the phonons to the left and to 
the right from the electron in Eq. J15II are interchanged 
(n m — ► n_ m ). If the two states \ti)k and \n) K are 
the same, Eq. ijlfil) should be replaced by the identity, 
\n + ) K = \n) K , while \u-)k = 0. 

The transformation given by Eqs. 11 (ill and Ijl7ll makes 
use of the point symmetry of the Hamiltonian (QJ , which 
is invariant under the space inversion. Because of this 
symmetry, the two subspaces ijlfil) and itTTjl are not cou- 
pled by the Hamiltonian for K = and K = it. 
Consequently, for K = and K = n the eigenstates 
have well defined parity properties: those belonging to 
the subspace iflfill are even, while those belonging to the 
subspace i!17|l are odd under the space inversion. Only 
the two polaron states of opposite parity can cross, as 
discussed in more details in Sec. IIVI 



III. NONADIABATIC LIMIT 

It is instructive to start the analysis of polaron prop- 
erties by inspecting the weak- and strong-coupling lim- 
its. In this way it is easier to understand the polaron 
formation in the regime for which neither of the param- 
eters in Eq. QJ can be treated as a perturbation. In 
Sec. IIII Al the ET results are compared with the weak- 
coupling perturbation theory (WCPT), while Sec. IIII Dl 
brings the discussion of the strong-coupling perturbation 
theory (SCPT). In the lowest order, the polaron disper- 
sion is described by nonadiabatic corrections for both 
these perturbative expansions. The intermediate regime 
of parameters (t ~ g ~ ujq), for which the polaron disper- 
sion involves a complex mixture of adiabatic and nona- 
diabatic effects, is investigated in Sec. IIVI In this latter 
case one has to rely on numerical results. 
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A. Weak coupling perturbation theory (WCPT) 

In the weak coupling limit the deformation field is 
treated as a perturbation to the free electron state. By 
applying the Rayleigh-Schrodinger perturbation theory 
the energy of the polaron states in the lowest band is, to 
the leading order in g, obtained as 

ik = 6 - 2 A/(wo-&)("o+4t-&) , (18) 

with £fc = 2t(l — cos(fc)). It is assumed in Eq. iflSjl 
that the polaron binding energy is small, i.e., g < too for 
t < ujq and g/wo < (t/wo) -1 ^ 4 for t > loqw^ This implies 
that the lattice distortion for any lattice site is less than 
the amplitude of the zero-point motion. Accordingly, one 
finds that, in the weak-coupling limit, the local polaronic 
correlations and the polaron dispersion are governed by 
nonadiabatic processes. In the absence of local adiabatic 
correlations, the spectrum below the phonon threshold 
consists of only one coherent band. 

Up to the second order in g, the phonon propagator © 
can be evaluated within the diagrammatic perturbation 
theory by considering the one-loop bubble diagram for 
the irreducible phonon self-energy. To the leading or- 
der in g, one obtains a phonon spectral function F q (ui) 
that, beside the phonon branch, exhibits two additional 
branches at ±£ q defined by the electron dispersioniSI For 
the latter case the spectral density is given by 



Notice that by considering only the first order Rayleigh- 
Schrodinger expansion of the polaron wave function |^ ) 
one obtains H q = 0. It is the second-order correction 
to l^o) that defines the lowest g 2 term in the expansion 
of H q . On the other hand, this term is taken into ac- 
count by the lowest-order diagrammatic expansion of the 
displacement-displacement correlation function^ 

Since the polaron dynamics in the weak-coupling 
regime is nonadiabatic, the static spatial correlations de- 
pend on the lattice mass. In particular, the g-dependence 
of Eqs. © and jf2T)l) is governed by the ratio t/ujQ. The 
electron-phonon correlation length d, can be defined as2& 

d 2 =J2 ?-2 (*o|<&c„ bl +r b n+r \^ ) I Ntot , (21) 

although other equivalent definitions are possible^^ 
With Ntot given by Eq. 112L one obtains d ~ yjt/uoo 
for t > loq and d ~ t/u>o for t < uiq. Irrespective of d, the 
second-order perturbative results in the weak-coupling 
limit are very close to the ET findings. This is illustrated 
by TableHl with d = 0.14 and d = 2.15 for t/u = 0.125 
and t/ujQ = 10, respectively. 





t = 0.125, g = 0.2 


t = 10, g = 1 


WCPT 


Ntot = 0.0272 


H = 0.0326 


Not = 0.08 


H = 0.156 


ET 


Not = 0.0275 


H = 0.0328 


Not = 0.083 


H = 0.159 



F q (uo = ±i q 



N (loo T t; q ) 2 



(19) 



The remaining spectral density belongs to the phonon- 
like branch at lo = loq, 

F q (u = cj ) = 1 + F q {u = -f,) - F q (uj = £,) . 
The last expression, derived here perturbatively, satisfies 



the boson sum- rule J1J. F q (w = £ q ) — F q (u 



< q ) de- 



fines the transfer of the phonon spectral weight from the 
phonon-like branch to = ujq to the electron-like branch 
lo = £ q . The additional spectral weight with respect to 
the free lattice is given by 2F q (ui = —£ q ). This is twice 
the mean number of phonons in the ground state © . 

When F q (oj) is known, H q (uj) can be determined 
by calculating the displacement-displacement correlation 
function ©. Up to the second order in g the spectral 
function H q (co) takes the form 



H q (L0) = 



9 



1 



Nu>i-Z 



S(lo ± £ q ) — — S(lo ± LOq 

LOO 



The integrated spectral density H q = J H q (io) duo is re- 
lated to the static correlation function 



H q = (* |M-,+&J&!j|*o>/2 = % , ( ,,\, s ■ (20) 



Table I: The second-order WCPT is compared to the ET 
method by means of Not and H = H q , where Not and 
H q are defined by Eqs. 11211 and 12011 . respectively. 

The lattice correlations in the frequency range \u>\ < uoo 
are investigated by the ET method in Fig. ^ f° r the 
t < loq and t > loq cases, using the same set of param- 
eters as in Table U For t/coo = 0.125 the coherent po- 
laron band, shown in Fig. QJi, is well-separated from the 
continuum of states above the phonon threshold (shaded 
area) for all momenta. The phonon spectra, shown in 
Fig. describe a polaronic lattice distortion which ex- 
hibits substantial quantum fluctuations, caused by the 
exchange of momentum between the electron and phonon 
subsystems. 

As it is well known for weak couplings^ the nature of 
the states in the coherent band below the phonon thresh- 
old changes as a function of momentum K if At > loq. 
This change corresponds to a crossover near K c , where 
K c denotes the momentum at which the phonon exci- 
tation energy is the same as the kinetic energy of the 
free electron, loq = £k c - For K < K c the momentum 
of the system is carried by the polaron (weakly dressed 
electron), while for K > K c it is carried by the phonon 
almost uncorrelated to the polaron in the ground K = 
state. In the latter case, the wave function of the sys- 
tem involves phonons at arbitrary large distances from 
the electron. Such long-range correlations, important for 
K > K c states in the lowest polaron band that are close 




ti/2 Q n ji/4 Q ji/2 



Figure 1: The lowest coherent polaron band and the lattice 
correlations in the weak-coupling regime obtained by the ET 
method. Panels a) and b) are for t = 0.125 and g = 0.2, 
panels c) and d) are for t = 10, g = 1, The shaded area in 
panels a) and c) corresponds to the continuum of states above 
the phonon threshold that characterizes the exact spectrum 
of the system. For the ET states that because of a small 
numerical error are located above the phonon threshold, a 
part (K > 7r/3) of panel d) is shaded. 



to the phonon threshold, are not described accurately by 
the ET approach. Although the exact lowest band re- 
mains below the phonon threshold for all K, the energy 
of the ET solutions for t/w Q = 10 in Fig. QJ for K > tt/3 
(K c ~ 7r/10) is slightly greater than the phonon thresh- 
old energy. In order to obtain better results one should 
use a set of basis states 11511 that includes phonons at 
larger distances from the electron than currently consid- 
ered. 

The lattice properties are analyzed in Fig. \Qi f° r 
t/ojQ = 10 by means of ET states satisfying EjP < 
Eq + u>q, whereas the shaded area corresponds to the part 
of the ET spectrum obtained above the phonon thresh- 
old (K > tt/3). The phonon spectral weight associated 
with the lowest band diverges in Eq. ijlQjl for Q — > K c 
as (ujq — £q) 2 - This singular behavior, replaced by an 
anticrossing within the degenerate perturbation theory, 
explains the rapid increase of F+q in Fig. QJi. On the 
contrary, the contributions related to the lattice distor- 
tion F_q vanish as K c is approached, with K~ x of the 
order of ground-state electron-phonon correlation length 
(for 4t > ujq). As shown in Fig. [Hi, with the exception of 
the long-wave limit Q«fl, the phonon softening effects 
are particularly intense, F+Jj F_q, in the narrow 
interval Q < K c . This weak-coupling phonon softening, 
associated with the long-range nonadiabatic correlations, 
should be distinguished from the softening of the local 
phonon modes. The latter, obtained for stronger cou- 
plings, is due to the short-range adiabatic correlations in 
the first place. 



B. Strong-coupling perturbation theory (SCPT) 

The small polarons are characterized by a deformation 
field localized mostly at one lattice site. For small bare 
electron bandwidths t < ujq the electron-phonon corre- 
lation length stays short for all couplings, whereas for 
t > ivo it is short when the coupling is sufficiently strong 
e p ^> t (A ;» 1). The two regimes, t < uiq and t > ujq, are 
commonly referred to as the nonadiabatic and the adi- 
abatic regime, respectively, although such identification 
oversimplifies somewhat the polaron problem. In both 
cases, the small polarons are frequently discussed using 
the SCPT. Within the SCPT, it is assumed that the elec- 
tron dispersion e(k) can be treated as a small perturba- 
tion of the polaron solution obtained in the atomic limit. 
The analytical expressions have been derived to the first 
and second order in t for the infinite lattice^ while the 
third order contributions are analyzed in Ref. US, but 
only for the two-site problem. From the low-order ex- 
pansion in t one is led to believe that the small polarons 
are necessarily heavy quasiparticles. Indeed, in the adi- 
abatic case one obtains small polarons strongly pinned 
by the discrete lattice, while in the nonadiabatic case 
the polaron is heavy because the bare electron hopping 
is small. It is argued here that, for t < ujq, the low 
order SCPT clearly demonstrates the dual character of 
small-polaron correlations: the nonadiabatic dispersion 
(translation) and the adiabatic softening associated with 
the local, adiabatic polaron deformation field. The mix- 
ing of nonadiabatic and adiabatic effects becomes par- 
ticularly important for i > wo, significantly reducing the 
small-polaron effective mass. However, the appearance of 
the light small polarons is beyond the reach of the low- 
order expansion in t (or g in the context of the WCPT) 
and therefore analyzed numerically by the ET method in 
Sec. |E| 

In the absence of hybridization, t = 0, there is no cor- 
relation in the lattice dynamics between different sites. 
At the site with the electron, the deformation field is 
described by a displaced harmonic oscillator, 

\<po) = c} e" 5 ' |0) , S 3 = -g/uo (fit - b,) , 

with the (zero-order) binding energy £o = ~£p- To the 
first order in t the polaron wave function takes the trans- 
lationally invariant form, 

\i Pl (k)) = l/VNj2 eikjc U~ Si \°) > ( 22 ) 

3 

corresponding to the well-known Lang-Firsov small-pola- 
ron band, £\{k) = —e p + e~ 1 e{k), where 7 = £ p /u>o and 
e(k) = —2t cos(fc) is the free electron energy. The dis- 
persion Si(k) is determined purely by nonadiabatic pro- 
cesses: during the hopping to the nearest-neighbor sites 
the electron detaches itself from the deformation field. 
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The probability of this hopping is reduced with respect to 
the free electron case by the electron quasiparticle weight 



|^i(fc)l4)| 2 = e- 



(23) 



For t < uiq, the first-order SCPT yields a satisfactory 
description of the polaron spectrum. This may easily 
be checked from the weak-coupling side for t -C loq. 
From the second-order WCPT Q18JI one obtains then 
the same polaron energy as from the 7 < 1 SCPT, 
Ei(k) w — E p + (1 — 7) e(k). In addition to the lowest 
band, the SCPT predicts one excited polaron band be- 
low the phonon threshold^ According to the degenerate 
perturbation theory to the first order in t, the first excited 
band detaches from the bottom of the phonon continuum 
for 7 > 1,12 



.(1) 



(k) = -e p + uj -t e~ 7 (x/2 + 2/x) , 



with x = 37 - 1 + y/J&y- !)(7- 1) for k = °> while 
for k = tt one obtains x = 7. In Fig. EK the first-order 
SCPT low-frequency spectrum (dashed curves), consist- 
ing of two bands, when compared for t = loq/& to the 
exact ET spectrum (solid curves) shows only small devi- 
ations. The i = and i = 1 bands, shifted by the ground 
state energy E , are plotted in terms of band boundaries 
(the K = and K = tt states). 

In the inset of Fig. Et the electron properties are in- 
vestigated by the ET method for g/uo = 1.52, while 
the lattice properties are investigated in Fig. Et>- The 
choice g/ujQ = 1.52 corresponds to the crossover regime 
between a nearly free-electron behavior, for weak cou- 
plings, and an exponentially-narrow-band behavior, for 
strong-couplings. The exact electron spectral density for 
the lowest i = band (the inset of Fig. EK) as a 
function of K differs slightly from the first-order SCPT 
value e -7 k 0.1. Similarly, the phonon spectral weight 
(Fig. El 3 ) for this band takes values close to Eq. (1111) . 
derived for the self-trapped small-polaron limit. In par- 
ticular, the three contributions of the i = band to the 
phonon spectra in Eq. Jl lll are obtained as J F 1 { 0) « 2.43, 

F! 0) w Ntot ~ 2.18, #(°) w 2.3. Because Fl 0) and N tot 
practically coincide, all the spectral weight Fq(u>) for 

lu < is assigned to the lowest band. Ff 1 — Fl ' w 0.25 
describes a small departure from the classical behavior 
(tl3ll of the local distortion in Eq. l(22l . This signalizes 
the interference between the polaron motion and the local 
correlations. For the first excited i = 1 band, the electron 
spectral function exhibits a strong independence. 
The phonon spectral weight associated to this excited 
band indicates that one internal phonon mode of the de- 
formation field is considerably affected by the adiabatic 
softening; F^ w 0.54 and f[ 1] w w 0. On the 

strong-coupling side A » 1 of Fig. Et the excited band 
asymptotically approaches the phonon threshold, which 
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Figure 2: The two lowest polaron bands below the phonon 
threshold are plotted in the upper panel (a) as functions of g 
for t = 0.125. The first-order SCPT results are given by the 
dashed curves, while the solid curves are the ET results. The 
polaron dispersion is taken with respect to the polaron ground 
state energy. In the inset of the upper panel the electron 
spectral function A% is shown for g = 1.52, while the lower 
panel (b) shows the phonon spectral functions, F_£ q, F_^q 
and H 



makes the softening effects particularly small. Such be- 
havior is obtained perturbatively in Ref. |^. 

For t > uj the SCPT theory converges slowly. The 
exception is the A » 1 regime, in which the adiabatic 
contribution to the polaron dynamics is reduced by the 
very small size of the polaron. Up to the second order in 
t the polaron energy is^ 



s 2 {k) = -e p + e- 2 M0) t/uo [Ei (2 7 ) - 7o - In (2 7 )] 
+ e^s(k) + e- 2 M2fc) t/uo [Ei (7) ~ 7o - In (7)] , 

where Ei( 7 ) is the exponential integral, and 70 is Euler's 
constant. In the limit 7 > 1, the fc-dependent terms 
decay exponentially (the first-order nearest-neighbor and 
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the second-order next-nearest-neighbor hopping), while 
the polaron binding energy becomes strictly adiabatic, 

e 2 « -e p (1 + 1/A 2 ) . (24) 

The second-order energy gain in Eq. 1241) is related to 
the spreading of the adiabatic deformation field. Namely, 
unlike in Eq. J22I) . where the deformation field develops 
at the electron site only, the second order wave function 

n,m j 

involves the electron-phonon correlations at three sites. 
Apparently, this spreading should also affect the polaron 
dispersion in higher orders. Indeed, for pinned small po- 
larons it was found by the ET methodic that, in the 
regime t ~ g, the lowest band dispersion is described 
accurately by e -7 a ^e(k). Here a(X) < 1, which is a 
function of A only, represents the adiabatic correction to 
the nonadiabatic first-order hopping. 



IV. ADIABATIC ASPECTS OF SMALL 
POLARONS 

The polaron self-trapping marks the change from the 
free-particle- like to the exponentially-large effective mass 
behavior of the polaron. As shown in Fig. |2ji, for t < ujq 
the polaron dispersion is reasonably well given by the 
nonadiabatic first-order SCPT for all couplings, including 
the self-trapping crossover at g ^ luq- 

For t > wo, the physical interpretation of the polaron 
self-trapping constitutes a difficult problem because of 
the mixing of the nonadiabatic and adiabatic effects. It 
was first suggested by Eagles^ that for t > ujq the self- 
trapping may be described by the anticrossing of two 
qualitatively different, nonadiabatic and adiabatic, po- 
laron states. The first is a large light polaron state, sta- 
ble on the weak-coupling side. Increasing the coupling, 
the energy of such large polaron state becomes compa- 
rable to the energy of the small heavy polaron, pinned 
to the discrete lattice. The latter is a stable solution 
on the strong-coupling side. The sharp, but continuous 
self-trapping crossover is obtained in terms of the ma- 
trix elements involved in the hybridization of the large 
(nonadiabatic) and the small (adiabatic) state. 

One might expect (for short-range electron-phonon in- 
teractions) that such an anticrossing scenario for t > ujq 
explains correctly the polaron self-trapping for two or 
more dimensional systems, in which the adiabatic po- 
laron is stable only when it is small j^ffl exhibiting 
strong pinning effects. In such circumstances the polaron 
self-trapping coincides with the nonadiabatic-adiabatic 
crossover. In ID, however, the situation is different. 
In particular, the regime of small adiabatic polarons is 
reached as g increases by two crossovers in the nature 
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0.5 1 1.5 2 A, 



Figure 3: The 8 lowest normal modes of the polaron defor- 
mation, obtained by the harmonic adiabatic approximation, 
as functions of A = e p /t. The pinning and the breather mode 
are denoted by top and ujb, respectively. 




Figure 4: For t = 5 the ET polaron spectrum below the 
phonon threshold, consisting of the lowest and three excited 
bands, is shown for the crossover regime as a function of g. 
In the inset the exact ET ground-state binding energy A po i 
is compared to the second-order WCPT and SCPT results. 



of the ground state, neither of which involves anticross- 
ing among the eigenstates of the Hamiltonian. That 
is, in the limit t ^> u the large nonadiabatic weak- 
coupling polaron crosses into a large adiabatic polaron for 
y/ujo/t r~j AiSI The polaron self-trapping involves a differ- 
ent crossover that occurs for A ~ 1. For \fu)o/t < A < 1 
the adiabatic polaron of the width d ~ A -1 behaves prac- 
tically as a free particle, because the pinning potential 
is exponentially smallf2& Within the harmonic adiabatic 
approximation^ the crossover from the free-particle- like 
to the pinned states is reflected in the disappearance of 
the restoring force near Awl for the lowest odd pin- 
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ning mode, as seen in Fig. In summary, the polaron 
self-trapping in ID can be caused by two different mech- 
anisms. The nonadiabatic one dominates in the t -C ojq 
limit. It can be described as the simple electron dressing. 
The other polaron self-trapping mechanism, dominating 
in the t ujq limit, is entirely adiabatic. 

The particularity of the moderate values of adiabatic 
ratio t/cu investigated below is that in the crossover 
regime between weak and strong couplings a mixture 
of two self-trapping mechanisms, nonadiabatic and adi- 
abatic, occurs. In particular, one finds significant nona- 
diabatic contributions to the polaron dispersion accom- 
panied with strong adiabatic softening of the two lo- 
cal modes of the deformation field, the latter indicating 
the weakening of the pinning potential that character- 
izes the adiabatic polaron motion. The ET spectrum 
in the crossover regime for t/ujQ = 5 is shown in Fig. 0] 
as a function of g. Only the states below the phonon 
threshold are displayed, shifted by the ground state en- 
ergy Eq. The lowest and the three excited polaron bands 
are represented each by the 22 states with different mo- 
menta, K = 0, 7r,n x 0.15, with n < 20. The spectrum 
in Fig. 31 reveals two important effects of the polaron 
motion. First, when the bandwidths increase with de- 
creasing g, there is a strong hybridization of the first 
and second excited bands. Second, when g is decreased 
further all the excited bands shift towards the phonon 
continuum. This latter effect indicates that the nonadi- 
abatic dynamics prevails in this limit. The fast electron 
gains more energy by detaching itself from the deforma- 
tion field, than by waiting to move adiabatically with the 
lattice distortion to the neighboring site. 



A. Strong pinning 

The very narrow polaron bands on the right side of 
Fig. 01 describe self-trapped small-polaron states. The 
bandwidths W^' are exponentially reduced by the strong 
pinning potential due to discreteness of the lattice. The 
positions of the excited bands in the spectrum are well 
predicted by the harmonic adiabatic approximation^ 
shown in Fig. 03 Because the polarons are small the 
adiabatic softening is notable only for short-ranged cor- 
relations, i.e., for the even (breather) and odd (pinning) 
modes involving three lattice sites. 

The self-trapped polaron represented in terms of Wan- 
nier states, 



E' 

K 



(25) 



stays in the quantum state \ifn ) centered at the lattice 
site n over a long period of time, l/W^ ^> l/uj a . In such 
circumstances the spatial (Q-dependent) lattice correla- 
tions are almost unaffected by the polaron dispersion, 
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Figure 5: The ET phonon spectral weight is shown as a func- 
tion of Q for A = 4 and s p /ojo = 10. Ntot is the mean number 
of phonons in the ground state 11211 . The thin dotted curves 
are given by Eqs. (23 and (2HJ- 



i.e., 



7 (i) 
+,Q 



(26) 



As obtained in Eq. ltTT)l . in the limit A> 1 the lattice is 
distorted significantly at one lattice site only, which in the 
reciprocal space results with the roughly Q-independent 
spectral weight of the lowest polaron band. However, 
this behavior changes by decreasing A as the distortion 
spreads to neighboring sites. For A = 4 and e p /uo = 10 
the phonon spectra are plotted in Fig. The parame- 
ters used correspond to the self-trapped polaron regime, 
in which F^q, ^ + qi an d Hq^ are practically the same 
(classical distortion). Hence, the decrease of the i = 
spectral weight toward the end of the Brillouin zone, 
which is seen in Fig. is related solely to the spreading 
of the local distortion. This spreading is characterized by 
the adiabatic correlation length d ~ 1/A as long as the 
distortion behaves classically. Increasing the range of the 
adiabatic correlations (at smaller A) only makes the in- 
dependent feature of the i = spectral weight observed 
in Fig. Elmore pronounced. 

In the A ^> 1 limit, the two lowest internal soft- 
phonon modes of the small-polaron deformation field 
involve three lattice sites, with a displacement pattern 
(1,-2,1) for the even (breather), and (-1,0, 1) for the 
odd (pinning) modei^l The corresponding phonon oper- 
ators with the appropriate normalization are given by 



, (i) _ b n -i — 2b n + b n+ i 



Q_n+1 ~ b n 



If the adiabatic fluctuations of the electron in Eq. (2fijl 
are neglected, the contribution of which is of the order 



10 



wo/£p, the phonon spectral density for the breather mode 
is obtained as 

* K°lWl°>l a = 3^ (l-cos(Q)) 2 , (27) 
and for the pinning mode as 

F[%*\(0\b^ Q \Q)\ 2 = ^ sin 2 (g). (28) 

Here, the lattice distortion is not considered explicitly as 
it does not contribute to F+'q . The expressions ffijl 
and (l28l explain the spatial correlations obtained by the 
ET method in Fig. for the two lowest i = 1,2 ex- 
cited bands of the self-trapped small polaron. Namely, 
the thin dotted curves, representing Eqs. ll2"7j) and p8|. 
compare well to the ET results for F^q (breather) and 

F^ Q (pinning). 

In contrast to the i = 0, 1, 3 bands in Fig. 0] the spec- 
trum of the second i = 2 excited band is inverted. The 
bottom and the top of this band are associated to the 
K = 7r and the K = state, respectively. Such dis- 
persion is a consequence of the parity properties of the 
local adiabatic correlations. This is easily understood by 
analyzing the self-trapped small-polaron limit, in which 
the hybridization between different bands is negligible. 
In this case, the Wannier states lj25)l inherit the sym- 
metry of the internal phonon modes of the deformation 
field. The pinning (i = 2) mode is odd under the space 
inversion, which results in a positive nearest-neighbor in- 
traband hopping. 

B. Crossover regime 

The bandwidths of the self-trapped polaron bands in- 
crease by decreasing the coupling due to the weakening 
of the pinning potential. The pinning potential is more 
effective for the polaron states at the bottom of the spec- 
trum, explaining why the bandwidths of the two lowest 
excited bands on the right side of Fig. 31 are significantly 
larger than the bandwidth of the lowest band. In the 
crossover regime, corresponding to the central part of 
Fig. Q] the bandwidths and the positions of the bands in 
the spectrum are characterized by the same time scale. 
This means that the local correlations, dominated in the 
self-trapped polaron limit by the harmonic adiabatic dy- 
namics, in the crossover regime are strongly affected by 
anharmonic and nonadiabatic contributions introduced 
by the small-polaron hopping. 

In the central part of Fig. Q] the polaron spectrum 
is characterized by the strong hybridization of the first 
(i = 1) and the second (i = 2) excited bands, while 
the hybridization between states of the second (i = 2) 
and third (i = 3) excited bands is weak. The i = 1 
and i — 2 bands touch for g ~ 3.32. This band touch- 
ing corresponds to the symmetry allowed crossing of the 
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Figure 6: The phonon spectral weight q, as a function of 
Q, for the first i = 1 (dashed curves) and the second i = 2 
(solid curves) excited band. The results are obtained for the 
set of couplings that corresponds to the narrow interval in 
which the two i = 1 and i = 2 excited bands touch (2.82 < 
g„ < 2.89, with g n +i = g n + 0.01, t = 2.5). The arrows point 
to the direction of the increasing coupling. The thick dashed 
and solid curves are plotted for the self-trapped small-polaron 
limit in Eqs. J23 and 

two K = 7r states of opposite parity. The origin of the 
strong hybridization of the two lowest excited bands lies 
with the small-polaron hopping. While the motion of 
the large polaron to the neighboring site corresponds to 
a small change of the distortion, the small polaron moves 
by creating the lattice distortion at one site and destroy- 
ing it at the other. The displacement pattern involved in 
the nearest-neighbor small-polaron hopping to the right 
is approximately given by (0,-1,1). It is easy to see 
that such a pattern is obtained by a linear combination 
of the breather (1,-2,1) and pinning (—1,0,1) displace- 
ment patterns. Thus, for small polarons the breather 
and the pinning mode are, due to discreteness of the 
lattice, effectively coupled by the polaron motion. Fig- 
ure H3 demonstrates the correlation in the behavior of the 
phonon spectral weight for the i = 1 and i = 2 bands 
in the touching regime. Figure El shows that, as the hy- 
bridization with the increasing coupling becomes weaker, 
there is a rapid transfer of the breather-like phonon spec- 
tral weight near the end of the Brillouin zone from the 
i = 2 band (solid curves) to the i = 1 band (dashed 
curves) . 

The polaron density of states (DOS), 

S(l>) = l/N Sk(uj) , S K (u>) = J2 ~ E k) . 

K i 

the phonon density of states (PDOS), F(u>) = 
l/N J2qFq(w), and the electron density of states 
(EDOS), A(oj) = l/N Y<k A k{u), are investigated in 
Fig. [3 for the crossover regime (t/u>Q = 5 as in Fig. 
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(@J). The left column of panels presents the coher- 
ent low-frequency features at the strong-coupling side of 
the crossover regime, g/coo = 3.35, whereas the weak- 
coupling side of the crossover regime is analyzed in the 
right column, g/uo = 2.9. For the former case the po- 
laron effective mass is large, m po iJm e i ~ 100, while for 
the latter case it is within the order of magnitude of the 
electron effective mass, m po i/m e i w 5. 

The DOS in Fig. contains three distinct features, 
each corresponding to one of the three lowest polaron 
bands in Fig. Q] (the contribution of the third excited 
band located in Fig. 21 very close to the phonon thresh- 
old is not shown in Fig. 0. Three distinct contribu- 
tions may also be observed in the PDOS of Fig. Et>- 
The low frequency feature, u) ~ in Fig. Et> is 

given by the narrow, <C ujq, lowest band, for which 
the phonon spectral weight behaves approximately classi- 
cally, F(u>) w F(—cj). Because of the weak Q-dependence 
of the phonon spectral weight (small-polaron distortion) 
the PDOS and DOS for the lowest band have similar 
shape. The considerable flattening of the lowest polaron 
band toward the end of the Brillouin zone is the reason 
for the more pronounced singular behavior of the DOS 
(and the PDOS) at the upper rather than at the lower 
band edge. 

The two lowest excited bands, related to the pin- 
ning and breather modes of the deformation field, al- 
most touch in Fig. EH- Their bandwidths are signif- 
icantly larger than the bandwidth of the lowest band, 
1^(1) -|- 1^(2) ~ (jj /2. The phonon spectral weight as- 
signed to the excited bands, accruing only for u> > 0, 
exhibits a strong Q-dependence. This explains the sub- 
stantial difference in the shape of the DOS in Fig. UK and 
the PDOS in Fig. 03 for frequencies corresponding to the 
1 = 1 and i = 2 bands. The hybridization of the breather 
and pinning modes due to the polaron hopping can be 
observed in the partial transfer of the phonon spectral 
weight from the i = 1 to the i = 2 band. 

For g/ujQ = 3.35 the electron properties are shown in 
Fig. Et- Although the polaron is small and heavy, the 
EDOS of the lowest band differs essentially from the sim- 
ple momentum-independent Lang-Firsov formula II23H . 
The most significant contributions to the EDOS involve 
the i = and i = 1 states near the bottom of the bands. 
In fact, one finds that, in the crossover regime, the elec- 
tron spectral weight for the lowest band is weakly de- 
pendent on the momentum only for small bare electron 
bandwidths (Fig. EJ. The contribution of the i = 2 band 
to EDOS in Fig. Et is very small because of its par- 
ity properties. While the free electron state belongs to 
the subspace 11 (ill , the main contribution to the polaron 
states of the i = 2 band is given by the states belonging 
to the subspace ifTTjl . In particular, the K = and K = it 
states of the i = 2 band are odd under space inversion, 
with no overlap whatsoever with the free electron state 
of even parity^ 

For g/oJo = 2.9 the stable polarons are small and light, 
exhibiting particularly interesting properties. Figure 01 
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Figure 7: The polaron (DOS), phonon (PDOS), and electron 
(EDOS) density of states for frequencies below the phonon 
threshold. The first and the second column of panels are 
given for g = 3.35 and g = 2.9, respectively, with t = 5 in 
both cases. The amplitude of the PDOS exceeds the range of 
values shown in Fig. 05. 



shows the polaron spectrum below the phonon thresh- 
old consisting of two coherent bands. The bandwidth of 
the lowest band is larger than that of the excited band. 
The PDOS is shown in Fig. Et- For the lowest band 
F(lu) shows notable asymmetry in ui, which means that 
the local lattice distortion begins to exhibit significant 
quantum low-frequency fluctuations associated with the 
polaron motion. Such quantum behavior is an indication 
of nonadiabatic effects. The adiabatic effects are mani- 
fested as a softening of the internal phonon mode of the 
deformation field. This softening is represented in Fig. Et 
by the excited band contribution, found only for uj > 0. 
The EDOS in Fig. is, in a manner similar to that seen 
in Fig. Et; characterized by two pronounced singularities 
at the bottom of the polaron bands. For g/too = 2.9 (as 
opposed to g/ujQ = 3.35) the polaron states are lighter 
for the lowest band than they are for the excited band. 
Accordingly, one finds that the lowest band contributes 
more to the electron spectral weight than the excited one 
in Fig. 0. Such a dispersion can be is explained by nona- 
diabatic dynamics that is energetically more favorable for 
the polaron motion than the adiabatic one. 





A po ; 


d 


Ntot 


m e i/m po i 


g = 3.35 


3.69 


0.45 


7.56 


0.01 


5 = 2.9 


2.25 


0.93 


2.25 


0.196 



Table II: Comparison of small-polaron properties for g = 3.35 
and g = 2.9 (t = 5). 

Table |TT| brings some additional insights upon the 
choice of parameters used in Fig. A po ; is the bind- 
ing energy, the electron-phonon correlation length d is 
defined by Eq. <EU, N tot is given by Eq. lfT2l . and 
m e i/m po i determines the polaron effective mass. By go- 
ing from g/uj = 3.35 to g/uJo = 2.9 the static polaron 
properties (A po ;, d, N to t) change rapidly, but far more 
dramatic is the change in the dispersion. For g/ujQ = 2.9 
the polaron effective mass is within the order of magni- 
tude of the electron effective mass, the latter being char- 
acterized by the large bare electron bandwidth, t/u>o = 5. 
It is emphasized that for such a large ratio of t/u>o, the 
motion of the polaron in the crossover regime involves 
substantial nonadiabatic and adiabatic effects. In the 
nonadiabatic case, the hopping occurs by a dephasing 
between the electron and the lattice distortion. The adi- 
abatic dynamics contributes to the small-polaron disper- 
sion through the softening of the internal phonon mode, 
which for small polarons, is obtained as a mixture of the 
breather and pinning modes. Because of this adiabatic 
softening the effective hopping overlap between the po- 
laron states at neighboring sites is larger than for the 
unrenormalized lattice vibrations. 



V. SUMMARY AND CONCLUSIONS 

The Holstein polaron problem is analyzed here for the 
one-dimensional infinite lattice. New important insights 
into the polaron formation are obtained and the crossover 
between the light and heavy polaron states is clarified. 
Using the practically exact polaron states, obtained nu- 
merically by the ET method, the spatial and temporal de- 
pendence of electron and lattice correlations is examined 
for the coherent part of the spectrum below the phonon 
threshold. Special attention is paid to the properties of 
the phonon propagator and the related correlation func- 
tions. To our knowledge, this constitutes the first inves- 
tigation of such kind in terms of translationally invariant 
polaron states. 

In order to achieve a better understanding of nonadia- 
batic and adiabatic effects in the polaron dispersion and 
of the local correlations, the weak- and strong-coupling 
perturbative results and the results of the harmonic adi- 
abatic approximation are compared to those of the ET 
method. The polaron dispersion may involve both adia- 
batic and nonadiabatic processes. In the first case, the 
electron is described by a localized state that follows adi- 
abatically the distortion as it moves along the lattice. In 
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the second case, the electron hops away from the dis- 
tortion before self-localizing again. The local adiabatic 
correlations are responsible for the softening of the inter- 
nal phonon modes of the polaron. Each of the excited 
coherent polaron bands found below the phonon thresh- 
old corresponds to one of these soft phonon modes. 

The distinction between nonadiabatic and adiabatic 
contributions provides a qualitative description of the 
polaron crossover between the self-trapped and free-like- 
particle behaviors of the electron, in the intermediate 
regime t ~ ujq ~ g in which neither of the parameters 
can be treated as a small perturbation. For t ~ luo ~ g 
the crossover regime involves the small polaron states, 
for which the adiabatic motion along the lattice is gov- 
erned by a pinning potential due to the discrete lattice. 
It is argued that the adiabatic contributions to the effec- 
tive nearest-neighbor hopping integral of such polarons 
are enhanced by the softening of the two lowest (even 
and odd) internal modes. In return, the polaron mo- 
tion effectively couples these two modes, which results 
in a hybridization of the states in the first- and second- 
excited polaron bands below the phonon threshold. This 
hybridization is weak on the strong-coupling side of the 
crossover regime for which the excited bands are narrow 
and well-separated. Because the pinning potential affects 
the low-frequency states more, the lowest band is notably 
narrower than the first two excited bands. As the cou- 
pling is decreased, the bandwidths of all the bands in- 
crease progressively because the pinning potential weak- 
ens. In addition, the hybridization between the two low- 
est excited bands becomes strong. The shifting of all 
the excited bands towards the phonon threshold ujq on 
the weak-coupling side of the crossover regime indicates 
that the local adiabatic correlations are suppressed by 
the nonadiabatic contributions to the dispersion. In this 
regime of competing adiabatic and nonadiabatic effects, 
one finds that the small polarons are light. The binding 
energy A po ; is a few ujq, while the spectrum below the 
phonon threshold consists of two coherent polaron bands. 
The polaron dynamics becomes fully nonadiabatic in the 
weak-coupling limit, in which the coherent part of the 
spectrum below the phonon threshold exhibits only the 
lowest polaron band. The weak interaction between the 
electron and the lattice distortion is reflected in the small 
binding energy (A po/ < u ). 

The lattice correlations change remarkably among dif- 
ferent regimes of parameters. The behavior of the phonon 
spectra is particularly complex in the crossover regime, 
in which the nonadiabatic and adiabatic contributions 
mix. The crossover regime connects two essentially dif- 
ferent limiting behaviors obtained for strong and weak 
couplings, respectively. In the strong-coupling limit (self- 
trapped polarons) the phonon spectral weight at very low 
frequencies describes a classical lattice distortion, while 
for frequencies corresponding to the soft internal phonon 
modes, the spectral weight is typical for harmonic dy- 
namics. The small \uj\ contribution, defined by the dis- 
tortion in the ground state, corresponds to the additional 
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spectral weight that is found with respect to the free 
lattice. The phonon softening involves the transfer of 
the spectral weight from the unperturbed frequency ojo 
to lower frequencies. As the slow quantum motion of 
the self-trapped polaron only weakly affects the local dy- 
namics, the Q-dependence of the phonon spectral func- 
tions describes the local spatial correlations. On the con- 
trary, for weak couplings the polaronic distortions exhibit 
strong nonadiabatic quantum fluctuations due to the ex- 
change of momentum between the electron and lattice 
subsystems. In particular, as the lowest polaron band ap- 
proaches the phonon threshold, significant phonon soft- 
ening effects are associated with long-range nonadiabatic 



correlations. The phonon spectral weight which corre- 
sponds to this softening can be much larger than that 
related to the weak-coupling distortion, which is defined 
by a small number of phonons in the ground state. 
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